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Abstract. This paper studies a closed queueing network containing a hub 
(a state dependent queueing system with service depending on the number 
of units residing here) and k satellite stations, which are GI/M/l queueing 
systems. The number of units in the system, N, is assumed to be large. After 
service completion in the hub, a unit visits a satellite station j, 1 < i < fc, 
with probability pj, and, after the service completion there, returns to the hub. 
The parameters of service times in the satellite stations and in the hub are 
proportional to One of the satellite stations is assumed to be a bottleneck 
station, while others are non-bottleneck. The paper establishes the continuity 
of the queue-length processes in non-bottleneck satellite stations of the network 
when the service times in the hub are close in certain sense (exactly defined in 
the paper) to the exponential distribution. 
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1. Introduction 

In the present paper we study the continuity of large closed queueing networks 
with bottlenecks. The continuity of queueing systems is a known area in queueing 
theory, and the goal of this theory is to answer the question: how small change of 
input characteristics of the system affects its output characteristics. This question 
is of especial significance for queueing networks that models telecommunication 
networks. 

In the present paper the continuity results are established for closed queueing 
networks describing client /server telecommunication networks. The queueing net- 
works considered in the paper are of the following configuration. There is a large 
server (hub), which is a specific state dependent queueing system (the details are 
given later) and k satellite stations, which are single server GI/M/l/oo queueing 
systems. Being served in the hub, a unit is addressed to the jth satellite station, 

f " \ 

^ 1^ j 1^ k, with probability pj > ^ Pj = ^ ] 7 si,nd then, being served there, it 

\j=i J 

returns to the hub. 

There are many papers in the literature that study similar type of queueing sys- 
tems. In most of them the hub is a multiserver (or infinite server) queueing system. 
The earliest consideration seems to go back to a paper of Whitt [26^ . A Markovian 
queueing network with one bottleneck satellite station has been studied by Kogan 
and Liptser [18]. Extensions of the results 18 for different non-Markovian models 
have been made in a series of papers and a book by the author (see [T], [2], [3], [S] 
and [8j). Non-Markovian models of networks with a hub and one satellite station 
have also been considered by Krichagina, Liptser and Puhalskii [19] and Krichag- 
ina and Puhalskii [20] . For other studies of queueing networks with bottleneck see 
Berger, Bregman and Kogan [TU], Brown and Pollett [TT] and PoUett [21 and other 
papers. 

Aforementioned papers [1], [2], [3], [5] and [18], study non-stationary queue- 
length distributions in non-bottleneck satellite stations. Specifically, in [T], [3] and 
[l8] these distributions have been studied in the presence of one bottleneck satellite 
station, and in [2] the only cases of bottleneck leading to diffusion approximations 
have been investigated. One of the main assumption leading to substantial technical 
simplification in [I], [5], [3] and [TH] was that at the initial time moment all of N 
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units are located in the hub, where the time instant < = is the moment of service 
starts for aU of them. In the case of exponentiaUy distributed service times in the 
hub, because of the property of the lack of memory of an exponential distribution, 
the above assumption is not the loss of generality. The initial time moment t = is 
not necessarily the moment of service starts for the units presented at the hub, and 
the number of these units need not be N exactly; it can be a large random number 
Ao which is asymptotically equivalent to N, i.e. =^ converges in distribution to 1 
as N increases unboundedly. 

However, in the case where the service times in the hub are generally distributed, 
the above assumption is an essential restriction, because its change can lead to 
a serious change of the limiting non-stationary queue-length distributions. The 
models considered in [2] (see also [8], Chap. 2) are flexible to this change, since 
they model hub as a state dependent single-server queueing system. The behavior 
of the queue-length processes in that system has been studied analytically, and the 
analysis of that system in [2] and [H| shows that it is stable with respect to a slight 
change of the initial condition. The assumption on state dependence in [2] and [S] , 
Chap. 2, was quite general, and real bottleneck analysis for that network, except 
the cases where diffusion approximations are available, could not be provided. So, 
the most interesting cases of the problem of the bottleneck analysis in [2] or [8], 
Chap. 2 remained unsolved. 

In the present paper we study a more particular model than that in [2], i.e. we 
study the model with the same state-dependent service mechanism in the hub as 
in [2], but specific additional assumptions on probability distribution functions of 
service times in the hub that will be discussed later in the paper. One of these 
assumptions, saying that service times in the hub are "close" (in different senses 
given in the paper) to the exponential distribution, enables us to easily extend the 
results to various models, including the realistic cases where the hub is multiserver 
(or infinite-server) queueing system. 

Recall the assumptions that made in [2] in the description of the service mecha- 
nism in the hub. In j2j it was assumed as follows. 
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(i) if immediately before a service start of a unit the queue-length in the hub is 
equal to i^T < A^, then the probability distribution function is Fk{x) = Gk(Kx). 
It was assumed that = xdGK{x)- 

Hence, for the expectation of a service time we have 

/ xdFK{x) = / xdGK{Kx) = — xdGxix) 
Jo Jo ^ Jo 



K\ 



K 

(ii) As — !■ oo, the sequence of probability distribution functions Gn{x) was 
assumed to converge weakly to G{x) with G(0+) = 0, and \ = xdG{x). In 
addition, it was assumed in [2 that Gk{Kx) < Gk+i{Kx + x) for all a; > 0. 

In the present paper in addition to (i) we assume that 
(1.1) Gi{x) = G2{x) = ... = Gn{x) := G{x), 

and there exists the second moment 



/•oo 

/ a;^dG(a;) ~ r < oo 
Jo 



satisfying the condition: 

(1.2) r>A. 



Condition (|1.2I) will be used later in the proof of Theorem 14.11 

In this case, for the random variable (^k having the probability distribution 
function Fk{x) = G{Kx) we have 

poo poo -1 

(1.3) ECk = / xdFKix) = / xdG{Kx) = —. 

Jo Jo 

Under the above addition assumptions, earlier assumption (ii) of ^ becomes irrel- 
evant here. 

Following (HH]) and dO]) one can "guess" that departures of units from the hub 
when K increases indefinitely should be asymptotically insensitive to the type of 
the distribution G{x), i.e. the number of departures per unit of time should be the 
same as in the case when G{x) would be the exponential distribution. This guess 
however cannot be proved by the means of the earlier papers [IS and [1 and should 
involve more delicate methods (which are provided in the present paper). 
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In addition to the assumptions made, we consider two classes of probability 
distribution function G{x) := P{C < x}. The first one satisfies the condition 

(1.4) sup \G{x)-Gy{x)\ < e, 

x,y>0 

where e is a small positive number, and Gy{x) — P{(^ < x + y\Q > y}. 

For the second class of distributions along with (|I.4p it is assumed that G{x) 
belongs to one of the classes NBU (New Better than Used) or NWU (New Worse 
that Used). Recall that a probability distribution function S(a;) of a nonnegative 
random variable is said to belong to the class NBU if for all a: > and y > we 
have S(a; + y) < S(a;)S(2/), where S(a;) ~ 1 — S(x). If the opposite inequality holds, 
i.e. S(x + y) > S(a;)S(y), then S(a;) is said to belong to the class NWU. 

The assumption given by (jl.4l) involves many nice properties to the probability 
distribution function G{x). We review these properties now, but the details will be 
given then. 

First, according to (|1.4p the probability distribution G{x) is close to the expo- 
nential distribution in the uniform (Kolmogorov) metric in the sense that if e is 
small, then the distance in Kolmogorov's metric between G{x) and the exponential 
distribution with mean j is small as well (see later relations (|1.7p and (|1.10p V 

Second, according to assumption (i) the probability distribution of a current 
service time in the hub is not changed when a new arrival occurs from one of satellite 
stations to the hub during that service time. The change of service time distribution 
in the hub can be at the moments of service start only. Additional assumption (II. 4p 
enables us to consider more general schemes in which the change of probability 
distribution function of a service time in the hub can be taken into account also at 
the moments of arrival of customers to the hub from satellite stations. This means 
that the present model can be essentially extended, and the extension can include 
the models where the hub is a multiserver or infinite-server queueing system. 

The details for the explanation of the second property are given in (|1.5|) and 
()1.8|) below. Specifically, under assumption (|1.4p we have: 
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(1.5) sup \Gy,{x) ~ Gy^{x)\ 
x,yi>0 

< sup \G{x) ~ Gy^{x)\ + sup \G{x) ~ Gy^lx)] 

x.yi>Q x>0 

< sup \G{x) — Gy^{x)\ + sup \G{x) — Gy2{x)\ 

< 2e. 

Therefore 

(1.6) sup \Gy,ix) - Gy,{x)\ < 2e 

x,yi,y2>0 

as well. As well, according to the characterization theorem of Azlarov and Volodin 
[9], [6] we have: 

(1.7) sup|G(x) - (1 -e-^^)| < 2e. 

x>0 

If, in addition, the probability distribution function G{x) belongs either to the 
class NBU or to the class NWU, then, instead of (|1.6I) . for any y2 > we, respec- 
tively, have: 

(1.8) sup \Gy^{x) — Gy2{x)\ < sup max \G{x) — Gy.{x)\ 

< sup \G{x) - Gy,{x)\ 

x,yi>0 

< e. 

Hence, 

(1.9) sup \Gy^{x) — Gy2{x)\ < €. 
As well, instead of (|1.7p (see [3]) we have 

(1.10) sup|G(a;) - (1 -e-^^)l < e. 

x>0 

As in [2] , it is assumed in the paper that the service times in the satellite station 
J, = 1, 2, . . . , fc, are exponentially distributed with parameter N/ij, and /ij > Xpj 
for j — l,2,...,fc — 1 and fik < Apfe. This means that the fcth satellite station 
is assumed to be a bottleneck station, while the first fc — 1 satellite stations are 
non-bottleneck. The only this case is studied in the paper. 
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The problem of continuity of queues and networks is an old problem. It goes back 
to the papers of Kennedy jl7| and Whitt [25 , and nowadays there are many papers 
in this area. To mention only a few of them we refer Kalashnikov jl5l , Kalashnikov 
and Rachev [16], Rachev [22], Zolotarev [27], [28], Gordienko and Ruiz de Chavez 

The continuity results of the present paper are based on the estimates for the 
least positive root of the functional equation 

(1.11) z = G[^l-^lz), 

where G(s), s > 0, is the Laplace-Stieltjes transform of the probability distribu- 
tion function G{x) of a positive random variable having first two moments, and 
satisfying the condition /i xdG{x) > 1. It is well-known (see Takacs [H] or 
Gnedenko and Kovalenko [12]) that the least positive root of the aforementioned 
functional equation (|1.11|) belongs to the interval (0,1). Within the interval (0,1) 
Rolski [25 established the possible lower and upper bounds for the least positive 
root of equation (jl.lip . when the probability distribution function G{x) has two 
moments. Let 5 (g 1,02) be the class of probability distribution functions of pos- 
itive random variables having the first two moments Qi and 02- Let G'^^\x) and 
G^^''(a;) be two probability distribution functions from this class satisfying the ad- 
ditional condition sup^>o \G^^\x) ~ G^^\x)\ < e. Let G^^\s) and G'^^\s), s > 0, 
denote the Laplace-Stieltjes transforms of G^^^(a;) and G^^''(a;) respectively, and 
let 7(^-' and 7'^^ denote the corresponding least roots of the functional equations 
z — G'^^^fi — iiz) and z — G^'^^^l — ^z). The lower and upper bounds for 17^^^ — t'-^''! 
have been established in [7]. In the present paper we aimed to use these bounds 
to establish the desired continuity results for the nonstationary queue-length dis- 
tributions in the non-bottleneck stations of our network. 

The rest of the paper is organized as follows. In Section|2]the known results from 
jl8) for Markovian models in which the hub is infinite-server queueing system and 
there is only one satellite station, which is a bottleneck station, are recalled. In same 
Section 121 we prove that the same result holds true under the main assumption of 
this paper where the service times in the hub are state dependent and satisfy (jl.ll) . 
In Section |3l we extend our result to the case of our model where there are several 
satellite stations and one of them is bottleneck. In Section HI the main results of 
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this paper related to the continuity of non-stationary queue-length distributions in 
non-bottleneck satellite stations are obtained. Some technical auxiliary results are 
recalled in the Appendix. 

2. A LARGE CLOSED NETWORK CONTAINING ONLY ONE SATELLITE STATION, 



In this section we consider examples of simplest queueing networks in order to 
demonstrate that in the case of the model considered in the paper, the asymptotic 
behaviour of the queue-length process in the hub is the same as that in the case of 
the exponentially distributed service times in the hub, which was originally studied 



2.1. The case G{x) = 1 — e^"^^. In the case where G{x) = 1 — e"'*'^ we deal with 
the model considered in the paper of Kogan and Liptser |18j . 

Consider the Markovian network (see Kogan and Liptser jTE]) containing only 
one satellite station, which is assumed to be bottleneck. Let A{t) and D{t) denote 
arrival process to this satellite station and, respectively, departure process from this 
satellite station. Denoting the queue-length process in this satellite station by Q{t) 
one can write: 



The departure process D{t) is defined via the Poisson process as follows. Service 
times in the satellite station are assumed to be exponentially distributed with pa- 
rameter fiN. Therefore, denoting the Poisson process with parameter fj.N by S{t), 
for the departure process D{t) we have the equation: 



Jo 

where l{-4} denotes the indicator of the event A. 

To define the arrival process, consider the collection of independent Poisson 
processes Tri{t), i = 1,2, . . . , N. Then, 



WHICH IS BOTTLENECK 



in [18]. 



(2.1) 



Q{t) = A{t) - D{t). 



(2.2) 




(2.3) 
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For the normalized queue-length process (?Ar(t) ~ j^Q{t), from (|2.ip we have: 

(2.4) QNit) = j^A{t) - j^D{t). 

As -/V — ^ oo, both jfA(t) and j^D{t) converge a.s. to the hmits, and these lim- 
its are the same as the correspondent limits of j^A{t) and -^D{t) (see [TH] or 
[T] for further details), where A{t) and D{t) are the compensators in the Doob- 
Meyer semimartingale decompositions A{t) — MA{t) + A{t) and, correspondingly, 
D{t) — MD{t) + D{t) {MA{t) and Af £>(/;) denote local square integrable martingales 
corresponding the processes A{t) and D{t).) Therefore, 

(2.5) qN{t) - - ^5(<). 

The compensators A{t) and D{t) have the representations (for details see |18j): 

(2.6) A{t) =X f [N- Q{s)]ds, 
and 







(2.7) D{t) = uNt ~^iN [ l{g(s) = 0}ds, 

^0 



Therefore, from ([231), ([22]) and ([SJ]) 
(2.8) q{t) : = P- lim <zjv(i) 



'OO 



= A /" [1 - (7(s)]ds - fit 
Jo 

+H f l{Q{s) = 0}ds. 
Jo 



'0 

In the case < A, using the Skorokhod reflection principle and techniques of the 
theory of martingales one can show that l{Q{s) = 0} vanishes in probability for all 
s > as increases to infinity (for details of this see [H] or [T]), and from 
we arrive at 

(2.9) q{t)=x[ [1 ~ q{s)]ds - fit. 

Jo 

Taking into account that g(0) — 1, from (|2.9p we obtain: 

(2.10) (z(i) = (l-^)(l-c-^*). 
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2.2. The case of arbitrary distribution G{x). We will show below that the 
equation (|2.10p holds also in the case, when service times in the hub are gener- 
ally distributed and state dependent (for further details see Section [T]) but under 
the assumption that there is only one satellite station, which is assumed to be 
bottleneck. 

Let (5 be a small enough positive fixed real number so the inequalities appearing 
later in the form /i < A(l — c5) for specified constants c are assumed to be satisfied; 
the large parameter iV is assumed to increase to infinity. Let QN{t) denote the 
queue- length in the satellite station (which is a bottleneck station), and let Q^it) 
denote the queue-length in the hub {Qpf{t) — N ~ Qjv(i)). In a small interval [0, S) 
let us denote 

Qn[0,6) = sup Q^it) 
o<t<s 

and 

Q~[0,S)^ mi Q^{t). 

U<t<o 

Clearly, that Q^[0,5) — N, and, as cxd, the a.s. lower and upper limits of 

9Ar(^) = t^Qn{^) as — i> oo, must be between 1 and 1 — (A — /i)(5 for all t G [0, 5). 

The last result is explained as follows. Let aAr[0, 5) be the maximum arrival rate 
(that is maximum instantaneous arrival rate) in the interval [0, 5). Since Q%[Q., 5) — 
N, this instantaneous rate unIO, — AA^. This means that, as — > oo, then with 
probability approaching 1 the number of arrivals during the interval [0,6) is not 
greater than XN6, i.e. 

(2.11) pJiin,sup^^^M<A4 = l, 

where An[0,S) denotes the number of unit arrivals from the hub to satellite sta- 
tion. We use limsup^_^QQ rather than limjv-s-oo because interarrival times are not 
identically distributed, so existence of the limit is not evident. 

Let iSAr[0,(5) denote the number of service completions in the satellite station 
during the interval [0,6). Similarly to (|2.1ip . as A^ — > oo, we have 

which means that, as A^ — > oo, then with the probability approaching 1 the number 
of service completions during [0, S) becomes close to fiNS. 
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The fact given in (|2.12|) can be easily proved as follows. Take e — jj, where M 
is a sufficiently large number. As fi < X, the probability that a busy period of the 
GI/M/l/co queueing system is finite is j (i.e. it is less than 1). Therefore the value 
N can be chosen such that during the time interval [0, e) there is only a finite number 
of busy periods with probability 1, and during the time interval [e, S) there is only 
an unfinished busy period, and units are served continuously without delay with 
the rate fj,N per time unit. Therefore, as N oo, with probability approaching 
1 the number of service completions during the interval [0, S), Sn[0,6), is at least 
fj.N{6 - e) = As iV ^ oo, the sequence of values M = M{N) can be 

taken increasing to infinity and associated sequence e — e{N) vanishing. Therefore, 
as — >■ oo, the value iSAr[0,(5) becomes asymptotically equivalent to jiNS. So, we 
arrive at ([2T2|) . It also follows from ([2ll|) and (|2?T2l) that 

(2.13) Pr (liminf rl^^Ml > A[l - (A - fi)S]s\ = 1. 

Af— s-oo N J 

From (P1T|) and ((^1^ wc obtain 

Pr (l - hmsup - ^) < liminf ,,(t) 



N N I ~ N^oo 



(2.14) < limsup9^(t) < 1 

N^oo 



= Pr <^ 1 - (A - < liminf g^(i) < limsupgjv(0 < 1 ^ = 1- 
Hence, according to (|2.14p 

(2.15) 1 - (A - fi)S '<q{t) <' 1 

for all t G [0,(5), where q{t) — i [liminf 7v-i.oo9Ar(i) + li™sup^^3^g^(t)]. Note, 
that as S vanishes, liminfAr^oo 97v(0 ^'^id limsupjv_j.o2 g^(i) approach one another 
(see rel. ([2TT|) and (f2TT3l) ). and from (|2T5l) we also have 

(2.16) linri^ = A-M. 

Sio S 

So, the value S can be chosen such small that q{6) is asymptotically close to 1 — 
(A — /i)(5. The existence of the similar limit is also obtained later in (|2.22p . (|2.24p 
and (12:26)) . 

Let us consider the second interval [S, 26). As A^ is large enough, for small S 
the number of units in the hub at time 6, Q]^{S), behaves as A^[l — (A — fJ.)S], i.e 



12 VYACHESLAV M. ABRAMOV 

we have approximately this number of busy servers in the hub. Therefore, at that 
time moment S the instantaneous rate of arrival of units in the satelhte station is 
asymptoticahy equivalent to AiV[l — (A — j.i)6] and the service rate fiN. Assuming 
the inequality /i < A[l — (A — /i)(5] satisfied, one can use the same arguments as 
above. First, qnI^, 2S) is asymptotically equivalent to AA^ [1 — (A — /i)(5]. As N ^ oo, 

(2.17) Prjlimsup^^^^^^^ < X5[l - {X - fi)S] \ = 1, 
and 

(2.18) p J linr ^^^Ml = ^5 

and at the moment 25 the queue- length in the satellite station, Qn{2S), is asymp- 
totically evaluated as ^[2{XS) - {XSf]N. It also follows from (I^TT)) and I^JE^ 
that 



(2.19) Pr|liminf:^^^I^>A<5 
From (|2Tfl) and (|2?T8| we obtain: 



1- V^[2(A5)-(A5)^] 



Pr 1 ^ limsup f — ^^-tt — ^ — ^tt — ^ I < liminf giv(i) 



(2.20) 



N N 



< lim sup <7jv (t) < 1 — lim sup 

A— >oo A— yoo 

Pril - ^^[2(A(5) - (A(5)2] < liminf g^(i) 

A A— >oo 



< Hmsup<7jv(*) < 1 ~ ('^ ~ /^)^ r ^ 1- 

A-i-oo 



Hence, according to (|2.20p 

(2.21) 1 - ^—^ [2{XS) - (A5)2] <'g(t) 1 - (A - ^)<5 

A 

for all t G [5^ 26). As S vanishes, similarly to f|2.16p we have 

Then, considering the third interval [26, 3(5), at the end of this interval the queue- 
length in the satellite station, QwiSS), is evaluated as ^[3{X5)-3{XS)^ + {XSf]N. 
In this case, similarly to (j2.2ip we arrive at the inequality 

(2.23) 1 - ^^[3(A(5) - 'S{XSf + [Xdf] '<g(^) <' 1 - ^^[2(A(5) - [XSf] 



CONTINUITY OF LARGE CLOSED QUEUEING NETWORKS WITH BOTTLENECKS 13 



for all t 6 [26,36). As 6 vanishes, similarly to (|2?T6| and (fZ22|) . 

(2.24) lim m)-m2SHms)-i ^ _ 



Thus, considering the sequence of intervals, at the end of the jth interval, for 
the queue- length in the satellite station, QN{j6), we obtain the expansion: 



X — fi 



J 



3 



{6X) 



i+i 



+{-iy {6xy] 

From this expansion we have 



(2.25) 




{6xy+^ + . . . + {~iy {5xy 



a.s a.s 

<g(i)<i 




for a\Ue[U-l)6,j6). 

Assume now that 6 vanishes. Then taking j — [|J, where [-J denotes the 
integer part of the number, and passing from the sequence of sums to the limiting 
expression, for q(t) we arrive at 

A — /i 



(2.26) 



m = 1- 

= 1 



A 

A — /i 
A 



l-lini(l-A(5)7 



(1 



-\t\ 



for all t e [0, oo). 

Summarizing all of this we have as follows. Consider the interval [0,oo), and 
build the system of partitions of this interval {H^}. Then for any point t G (0, oo) 
there is the integer value js (depending on 6) such that 6js < t < 5js +5. As 5 
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vanishes and, correspondingly, js increases to infinity, then in the point t we have 
the convergence q{t)^= Mhin^oo 9Ar(^), where q{t) satisfies (|2.26p . According to the 
construction, this a.s. convergence is uniform on any bounded subset of the interval 
[0,cx)). 

Hence, in view of q{t) = 1 —q{t), (|2.26p implies (|2.10p . and the desired relation- 
ship is proved. 



3. A GENERAL CLOSED NETWORK CONTAINING SEVERAL SATELLITE STATIONS, 

ONE OF WHICH IS BOTTLENECK 

Consider now a more general case of the model, where there are k satellite sta- 
tions numbered l,2,...,fc, and only the fcth satellite station is bottleneck, while 
the first k — 1 satellite stations are non-bottleneck. In this case as in P] we use 
the level crossing method. Note, that the approach of 2] based on straightforward 
generalization of the level crossing method related to Markovian queues was mis- 
taken. The error of ^2, was corrected in another paper [4 and then in book [8] , and 
here we follow by the improved method. Let us recall this method in the case of 
the standard GI/M/1 queueing system and then explain how the method can be 
adapted to this network with bottleneck. 

Consider the GI/M/1 queueing system, where an interarrival time has the prob- 
ability distribution function G{x) with the expectation j, a service time is exponen- 
tially distributed with parameter /i, and X < fi. For a busy period of this system, 
let f{i) denote the number of cases during that busy period where a customer at the 
moment of his/her arrival finds exactly i other customers in the system. Clearly, 
that /(O) — 1 with probability 1. Let U^i, ti^2, ■ ■ ■ , denote the time moments 

of these arrivals, and let s^^i, Si.2,. . . , ^^./(i) denote the service completions during 
the aforementioned busy period, at which there remain exactly i customers in the 
system. 

Consider the intervals 



(3.1) 
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and 

[ti+1,1, Si+l^l), [ii+l,2j •Si+l,2)i • ■ • , 

(3-2) Si+ij(i+i)) 

The intervals of (|3.2I) all are contained in the intervals of p.ip . Let us delete 
the intervals of p.2p from those of p.ip and merge the ends. Then we obtain 
a special type of branching process {/(«)}, which according to its construction is 
not a standard branching process, because the number of offspring generated by 
particles of different generations are not independent random variables. (Only they 
are independent in the case of the M/A//1 system.) However, the process satisfies 
the property: E/(i) — 93', i = 0, 1, . . ., which is similar to the respective property 
of a standard Galton- Watson branching process. The branching process {/(«)} is 
called GI/M/1 type branching process. The properties of this branching process 
are discussed in 4 and |8]. 

The value (p :— E/(l) is calculated as follows. We have the equation 

E/(l) ^^^Y.^^ r e-^-M:dG(x), 

where is the least positive root of the functional equation z = G{fi — fiz); G(s), 
s > 0, is the Laplace-Stieljies transform of the probability distribution function 
Gix). 

Let us now consider the closed queueing network with k satellite stations. The 
first k — 1 satellite stations are non-bottleneck, i.e. the relation fj,j > Xpj is satisfied, 
while the fcth satellite station is a bottleneck station, and then the relation /i^ < Xpk 
is satisfied (Recall that by pj , j — 1,2, ... , k, we denote the routing probabilities.) 
The input rate to the jth satellite station at time t is XpjMt, where the random 
variable Nt (A/j < N with probability 1) is the number of units it the hub at time 
t. 

By the level crossing method one can study the queue-length process in any 
satellite station j, j = 1, 2, . . . , fc — 1. Consider first the time interval [0, 5), where 
(5 is a fixed sufficiently small value. Then, similarly to the construction in Section [2] 
in the case of only one satellite station, for N increasing to infinity with probability 
approaching 1 the queue-length in the hub at time moment t — 5 is asymptotically 
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equal to N[l — {Xpk — IJ-k)S]. If there are several non-bottleneck satellite stations, 
then the queue-length in the hub at time moment S is asymptotically equal to the 
same aforementioned value Af[l — (Xpk — l^k)^] because the queue- lengths in non- 
bottleneck satellite stations all are finite with probability I, and their contribution 
is therefore negligible. More detailed arguments are as follows. 

In any satellite station j < k the length of a busy period is finite with probability 
1. Hence, as N increases unboundedly, the number of busy periods in the interval 
[0, 6) increases to infinity as well, and at time moment 6 the total number of units 
in all of satellite stations j < k with probability 1 is bounded, i.e. it is negligible 
compared to N. Therefore, the number of customers in the hub is asymptotically 
evaluated as A^[l — {Xpk — A'fc)(5] as well. 

For the jth satellite station, j < k, let us consider the last busy period that 
finished before time moment S. Let fj{i) denotes the number of cases during that 
busy period that a customer at the moment of his/her arrival find i other customers 
in the system, let ij,i.2,- • • ytj,i,fj{i) be the moments of these arrivals, and let 

Sj,i,2T ■ ■ ,Sj^ij.(^i) be the moments of service completions that there remain 
exactly i units in the satellite station. We have the intervals 

(3.3) [tjA,!, Sj^i^i), [ij,i,2, Sj,i^2), ■ • ■ , [^j,iJ{i)T^j,iJj{i)) 

and 

[tj,i+l,l, Sj [ij,i+l,2, Sj,i+1,2), • ■ • , 

(3.4) , Sj^i+ij. (i+i) ) 

which are similar to the intervals and p.2p considered before. (The only 

difference in the additional index j indicating the jth satellite station.) Delete the 
intervals p.4p from those of (|3.3|) and merge the ends. As N increases unboundedly, 
the process {fj{i)}i>o in a random sequence in i. Each of the random variables 
fj{i) converges a.s. (as N oo) to the limiting random variable, which is the 
number of offspring in the ith generation of the GI/AI/1 type branching process. 
So, we have the a.s. convergence to the GI/M/l type branching process in the 
sense that for all i this a.s. convergence holds. Let us find E/j(l). Similarly to 
the above case of the GljMjX queueing system for sufficiently large N we have the 
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equation 

1 - (1 -ft)^Lw[i-(Ap,-p,)5]j(/"j^ - ^^]Nz) 

(3.5) 



P3^lN[l-(\pk-t^k)S]\ y [N[l~(Xpk-fJ^k)S]i ) 
1 /"I \/^ { fj,jN~fijNz \ 



[I + 0(1) 



where Fn{s) denotes the Laplace-Stieltjes transform of the probabiUty distribution 
function F^ix) — G]y{Nx) = G{Nx). For details of the derivation of (|3.5p see the 
Appendix. As — oo, in hmit we have the equation 



(3.6) 



Considering now the interval [6, 26), in the endpoint 26, due to the arguments of 
Section [2] and the arguments above in this section, the queue-length is asymptoti- 
cally equal to N^^^^[2{Xpk6) - {\pk6)'^]. Therefore, similarly to ((5^ for the jth 
satellite station, j < k, we have the equation 



l-(l-p,)G 



l_ >-Pk^k ^2XpkSHXpkS)^ 



In an arbitrary interval [{i — V)6, id), in its endpoint i6 we correspondingly have 
the equation: 



(3.8) 



where 



PjG 



U{i,S) 



1-{1-P,)G 



U{i,S) 



+ (-1)' I ^ I (Apfc<5)' + ... 



-{-inxptsn 
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Taking i = [|J , in the limit as J — >■ we obtain 



(3.9) 



H__H± 

9(t) 



where q{t) = 1 - ^^^^{l - e-^^"*). 

Thus, denoting by ipj{t) the root of the functional equation p.9p . we arrive at 
the following statement. 

Theorem 3.1. The queue-length distribution in the non-bottleneck satellite station 
j < k is 

P{Q,{t) - = p,{t)[^,{t)Y-'[l - ^,{t)], 
(3.10) z = l,2,..., 



wh' 



ere 



M{t)pj 

Mi 



q{t) = l- ^P'' ^' (1-e-^P^*), 
Apfc 



and ipj {t) the root of the functional equation 

4. Continuity of the queue-length processes 
4.1. Formulation of the main results. 



Theorem 4.1. For small positive e assume that Condition (|1.4p is satisfied. Then, 
for ipj (t) in p.lOp the following bounds are true: 

Pj{t) - ei{t) < ifi.it) < p,it) + e2{t), 

where 

ei{t) = min (p,(t) - £,(t), -[1 - £j{t)] 

e,{t) = min |l + ^^^Ut) - 1] - " 
£j {t) is the least root of the equation 
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b,(t)- I- ^ ' ' ' > 



1 

Theorem 4.2. Under the assumptions made in Theorem \^.l\ assume additionally 
that the probability distribution function G{x) belongs either to the class NBU or 
to the class NWU. Then, for (pj{t) in (I3.10p the following bounds are true: 

Pjit) - esit) < ifijit) < p,it) + e4t), 

where 

eaW = min \ p,{t) - £,{t), -[1 - i,{t)] 
L Pj 

64(t) = min |i + ^^^im - 1] - P. ^[1 - m] 

and£j(t), aj{t) and bj{t) are as in Theorem \4-l\ 



Note, that the difference between ei{t) (or e2{t)) in Theorem 14. II and e3(t) (or 
e^it)) in Theorem 14.21 is only in one term in the min function. While ei{t) (or 
e2{t)) contains the term — £j{t)], the corresponding term in e3(t) (or 64 (t)) is 

4.2. Background derivations and the proof of the theorems. The results of 
continuity that used to prove Theorems 14.11 and 14.21 are based on a recent result 
obtained in [7 . Recall it. Let Gi{x) and G2(x) be two probability distribution 
functions belonging to the class Q{a,b) of probability distributions functions of 
positive random variables having the first and second moments a and b respectively, 
b > a^. Assume that 

(4.1) snp\Gi{x)~G2ix)\<K, 

where k < 1 — Let Gi{s) and G2(s) be the corresponding Laplace-Stieltjes 
transforms of Gi{x) and 6*2(2;) and let jci and 7G2 be the least positive roots of 
corresponding functional equations z = Gi(/i — fiz) and z = G2(/i — pz), where 
> i is some real number. Then \ — jq^ 1 < ~ ^)j where £ is the least root 
of the equation 
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On the other hand, according to the results of Rolski [23 , the guaranteed bounds 
for any of probabihty distribution functions Gi{x) and G2{x) (i.e. not necessarily 
satisfying (|4.1I) ) having the first two moments a and h respectively are given by 

a? 

(4.2) £<7G, <l + y(^-l), i = l,2. 
In our case, p.9p can be formally rewritten as 

z = H,{iij - ^ijz), 

where 

H(s) "-^(^) 

i-(i-p,)g(=^) 

is the Laplace-Stieltjes transform of a positive random variable having the proba- 
bility distribution Hj{x), 

oo 

G*'^{x) denotes z-fold convolution of the probability distribution function G{x) with 
itself, and^(i) = 1— ^^Xp^^ (1 — e^'^^''*). Our challenge now is to find the bounds for 
the least positive root of the functional equation in (j3.9p similar to those given in 
(|4.2[) and then, applying the aforementioned result of [7] , find the continuity bounds 
for non-stationary distributions in non-bottleneck stations. (The notation 7//^ is 
used here (along with the other notation ^j{t") in the formulations of Theorems 13.11 
14. II and [4. 2p for consistency with the notation such as and that introduced 
before.) 

The probability distribution function G[x) is assumed to have the expectation 
^ and the second moment r > -p- . Let us find the expectation and second moment 
of the probability distribution Hj{x). The best way for deriving these numerical 
characteristics is to use Wald's identities as follows. Let ^1, ^2, . . . be a sequence 
of independent and identically distributed random variables, let r be an integer 
random variable independent of the sequence ^1, ^2, ■ • ■ • Denote 5'r = ^1 + ^2 + 

+ Then, 

(4.3) ESr = E^Er, 

(4.4) Var(5^) = Var(a)ET -I- Var(r)(E^i)2. 
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In our case, the random variable r is a geometrically distributed random variable, 
Prjr — n\ ~ ■pj[\ — pj)"^^, n>\. Therefore, 

1 

Er= -, 



and 



Var(T) = + 



1 1 - 2pj 
2 ■ 



1 - Vj Pj 

The random variable has the probability distribution function G{q{t)x). Hence, 
and 

So, according to Wald's equations (|4.3p and (|4.4p we obtain: 

1 



(4.5) / xdHj{x) 
Jo 



OO 



(4.6)/ x^Hjix) = r 



1 / 1 \ / 1 l-2pj 



1 



AV[g(i)]2- 



Proof of Theorem \4-l\ It follows from (j4.5p and (|4.6p and the aforementioned 
result by Rolski [23] that satisfies the inequalities 

(4.7) ^,(i)<7H, <1 + M![^^.(i)_l], 

where (t) is the least root of the equation 

(4 8) ^ ^ ^-aj(t)^j+aj(t)/ijZ^ 



(4.9) a,(i) = 

Apjg(i) 

and 



(4.10) 



1 



+ 



x^pmt)?' 
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Assume that ()1.4p is satisfied. Then according to the characterization theorem 
of Azlarov and Volodin [9], [6] we have (|1.7|) . and then 
(4.11) 
G 



sup 

s>0 



Xq{t) 



q{t)J Xq{t) + s 



= sup 

s>0 . 


/"OO /" 

/ e-"="dG(g(t)a;) - / 
'o Jo 






< sup J 

s>0 J 


poo 

1 se-''' sup \G{q{t)x) - 

x>0 


-(1 


_e-^9(*)-)|da; 


= sup 1 

s>0 J 


poo 

1 se-''='sup|G(a;) - (1 

2;>0 


— e" 


~^^)|da: 




<2e 







< 2e. 

In turn, from (j4.1ip we obtain 

Xpjqit) 



sup 

s>0 



sup 

s>0 



Xpjq{t) + s 

oo 



(4.12) 



q{t)J \\q{t) + s 

f Xqjt) 
q{t)J \Xq{t)+s 



< 



-Pj) 



_ 2e 
Pj' 

Then, the results of [7^ enables us to conclude as follows. Let us consider the 
functional equation 

(4.13) ^-T^WT^ • 

Xpjq[t) + fij — fijZ 

The right-hand side of (|4.13p can be written as Il{p,j — p,jz), where n(s) is the 
Laplace-Stieltjes transform of exponential distribution with the parameter Xpjq{t) 
The least positive root of equation (|4.13p is 

Xpjq{t) 



Apparently, 
(4.14) 



Pi it) 



^At)<lH, < 1 + 
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where £j{t), aj{t) and bj{t) are defined in ()4.8p - (|4.10p . The lower and upper bounds 
in ()4.14|) are natural bounds obtained from the results by Rolski [23] . 
Along with (I4.14p we have also the following inequality 



(4.15) 



p,{t) - -[1 - tj{t)] < 7ff, < p,{t) - -[1 - e,{t)]. 
Pj Pj 



which holds true because of Condition p.2p . 

On the other hand, following the results in [7^ and (|4.12l) . 

\iH,~pAt)\<-[i-iM- 
Pj 

So, 
(4.16) 

Amalgamating (|4.14p and (|4.16p we arrive at the following bounds: 
(4.17) p,it) - eiit) < < Pj{t) + £2(0, 

where 

2e, 



em 



eiit) = min \ p,{t) - i,{t), -[1 - £,(1)] 

Pj 



and 

£2(0 = min <j 1 
The theorem is proved. 



Proof of Theorem \4-S\ The proof of this theorem is similar to that of Theorem 
14.11 The only difference that instead of (|1.7|) given by characterization theorem of 
Azlarov and Volodin [5] one should use (jl.lOp from |6j. In this case, instead of 
estimate (|4.1ip we arrive at 



sup 

s>0 



G 



and then instead of (|4.12p at 



sup 

s>0 



q{t)J \q{t) + s 



Xpjq{t) + s 



< 



Pj 



Then instead of (|4.16p we have 



p,{t) - -[1 - e,{t)] < < Pjit) - -[1 - iM^ 
Pj Pj 
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and instead of (|4.17p we obtain 
with 



e 



eaW = min \ pj{t) - £,{t), -[1 - £,{t)] 

Pj 



and 



12 e 



£4 W = min |l + ^^ij^Mt) - 1] - W, -[1 - m] 
The theorem is proved. 

Appendix: Derivation of (|3.5p 

Let Ti, T2, . . . be a sequence of service completions in the hub, and let QoiTi), 
Qo(j'2), ■ ■ ■ the queue-lengths in the hub in these time instants ti, r2, . . . . Then, the 
probability distribution function of the first interarrival time in the satellite node 
j can be represented as 

oo 
1=0 

where Fr„{x) = F^ix) = Gn{Nx), Fr,{x) = Y.u=i MQo(tO = u}Gu{ux), I > 1, 
the asterisk denotes convolution between the probability distribution functions. 

Denote the sequence of interarrival time distributions in the satellite node j by 
Hj^i{x), Hj^2{x),. . ■ . Let us derive the representation for Hj,2{x). Keeping in mind 
that the time instants ti, T2,. . . , r^, . . . can be the moments of the first arrival to the 
satellite node j with probabilities pj, pj{l — Pj),- . . , Pj{l — pjY~^,. . . respectively, 
we obtain: 

oo oo 
1 = 1 i=0 

and recurrently, 

oo oo oo 

(A.3) fc„_i=ii=o 

^ [-^'■fci+...+fc„_i * -^i-fci+.-.+fc^.i+i * ■ • • * -^'•fci+...+fe„_i+i(^) 
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Let US now use the assumption given by According to this assumption, we 

have Fi{x) < Fn{x) for all i < N. Hence, for the convolution given in the right- 
hand side of (A. 3) we have: 

where F^'^^{x) is {i + l)-fold convolution of the probability distribution function 
Fn{x) with itself. Therefore, 

oo 

(A4) Ha^) < Y^p,{l-p,yF»+\x) 

1=0 

for all w > 1 and a: > 0, and moreover formula (A. 4) remains correct when v is 
a positive integer random variable not smaller than 1, due to the formula for the 
total probability. 

Let us consider the intervals of the type of p.3p and p.4p for the node j. To this 
end, consider the procedure of deleting the intervals of p.4p from those p.Sp and 
merging the ends. Denote also and NJ' for maximum and minimum numbers 
of customers in the hub during the busy period in the jth satellite station where 
these intervals are considered. Clearly, that NJ < N, and Nj > N — 2_/i=i"'ii 
where denote the number of unit arrivals to the satellite node i from the hub. 
Let Vj denote the number of units served during the busy period in the satellite 
station j. The distribution of Vj generally depends on the state of the network at 
the moment of the busy period start, and it is rightly to write Vj{S). This state S 
includes the number of units in the hub, which is assumed to be large enough such 
that as — !> cxD it becomes asymptotically equivalent to A^. 

Apparently, for the expected number of units served during the busy period 
we have the inequality FiVj^S) < C, where C is the expected number of served 
units during a busy period of the GI/M/1 queueing system with independently 
and identically distributed interarrival times having the probability distribution 
^^QPj(l ~ pjYF^'^^{x) and exponentially distributed service time with parame- 
ter ^jN . (The mentioned probability distribution of a service time is given by the 
right-hand side of (A. 4).) For a non-bottleneck station j the value C is finite, since 
the expected interarrival time associated with the probability distribution function 
Y^'iLoPji^ PjY^N^^i^)^ which is equal to j^-j^, is assumed to be greater than 
— . Hence, the random sum X]i=i '^j is a finite random variable with the finite 



' ki + ...+k^_l ' ki+. . . + k^_i+l ' ki+ . . . + k^_i+i \ I 
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expectation. Therefore, as N increases to infinity, converges a.s. to f. This 
means that, as iV — > oo, the random variable i^j converges in distribution to the 
number of served units during a busy period of the GI /M/1 queueing system with 
independently and identically distributed interarrival times having the probability 
distribution function '^°^QPj{i—PjyF^^^{x) and exponentially distributed service 
times with parameter fijN. As well, let fj{S, n) denote the number of cases during 
the busy period i'j{S), when an arriving unit meets n units in the jth satellite sta- 
tion. As N oo, the process {/j(iS, n)} converges in distribution to the associated 
G//A//l-type branching process. Specifically, let fj^N be the least positive root of 
the functional equation 

oo 

= Y^Pi^^-P^y^N^t^N - ^lNz) 

i=0 



l-[l-p,)FM{tMjN - ^lJNz)' 



Then, 



lim YifAS.l)— lim zm 



N^ca 1 _ (1 _p^.)Fjv(AijiV - PjNz) 

Taking into account, that F/v(s) ~ Gn (■§■) and limjv^oo Gn{s) = G(s), from (A.S) 
we obtain 

lim E/,(5, 1) = = PMH^-i^m) , 

N^oo 1 - (1 - Pj)G{^lj^j - fljifij) 

where ipj = limjv^oo '^],n- 

Let us now derive (|3.5p . Let [0,5) be a small time interval, and let Z be the 
last busy period in this time interval in the satellite station j. As N becomes large 
enough, at the end of the interval [0, S) the number of units in the hub is proportional 
to [1 — {Xpk ~ fJ-k)^]N. The following arguments in the derivation of (|3.5p are the 
same as above. The only difference is that the probability distribution function 
F/v(x) is to be replaced with i^L[i-(Apfc-/^fc)(5]Afj and, consequently, its Laplace- 
Stieltjes transform Fn{s) is to be replaced with -F[[i-(Apfc-/^fc)5]A'j (s)- Therefore, 
instead of asymptotic relation (A.S), we obtain the following asymptotic relation 

E/,(5,l) = pAml-i^r>.-,Mii^^^N-^,,Nz) ^^^^^^^^^ 

1 - (1 -Pj)^'L^[l-(Apfc-M;=)5]j(Mi^ - l^jNz) 
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which coincides with the asymptotic formula in p.5|) . 
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